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This paper is concerned with all tests for continuous- variable entanglement that arise from linear combina- 
tions of second moments or variances of canonical coordinates, as they are commonly used in experiments to 
detect entanglement. All such tests for bi-partite and multi-partite entanglement correspond to hyperplanes in 
the set of second moments. It is shown that all optimal tests, those that are most robust against imperfections 
with respect to some figure of merit for a given state, can be constructed from solutions to semi-definite opti- 
mization problems. Moreover, we show that for each such test, referred to as entanglement witness based on 
second moments, there is a one-to-one correspondence between the witness and a stronger product criterion, 
which amounts to a non-linear witness, based on the same measurements. This generalizes the known product 
criteria. The presented tests are all applicable also to non-Gaussian states. To provide a service to the commu- 
nity, we present the documentation of two numerical routines, FullyWit and MultiWit, which have been 
made publicly available. 

PACS numbers: 03.67.-a, 03.65.Ud, 42.50.Dv 



I. INTRODUCTION 

The field of continuous-variable quantum information has 
seen a very substantial progress in recent years. This has in 
part been made possible - from the experimental side - by 
the availability of a number of sources of systems prepared 
in entangled states in the canonical coordinates. Notably, 
two-mode squeezed states of Ught close to minimal uncer- 
tainty have been prepared (7], as well as bright entangled light 
beams |2, 3]. The collective spin states of atomic ensembles 
have been brought into states that can be well-described in 
terms of continuous-variable entanglement 14*1, even allowing 
for a light-matter interface |5]. First instances of multi-mode, 
multi-partite entangled states have also been prepared already 

In many of these set-ups, the starting point of further ex- 
ploitation of entanglement is to see whether the envisioned 
bi-partite or multi-partite entanglement can be found in the 
prepared states. This is often done by making use of criteria 
of entanglement based on second moments, or uncertainties, 
of quantum states. Notably, in a two-mode set-up, the proba- 
bly most well-known criterion of this type is the following: if 
one finds that the state p of a two-mode system equipped with 
the canonical coordinates xi and pi of one mode and X2 and 
■p2 for the second mode fulfils 

{{u-{u)p)'), + {{v~{v)p)')p<l (1) 

where u = {xi +X2)/ \/2 and v = {pi —P2)/ v^, then one can 
assert that the state must have been entangled Such a crite- 
rion is tremendously helpful: firstly, it gives a clearcut test for 
deciding whether a state is entangled or in a subset where one 
cannot assert whether it was separable or entangled. Secondly, 
one only has to measure certain fixed combinations of sec- 
ond moments of the original canonical coordinates. This is, 
yet, only one specific test, detecting the entanglement in some 
states, not detecting it in others. Similar tests have been pro- 
posed also to detect the entanglement of certain multi-party 
entangled states 
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FIG. 1: (a) Schematic representation of the optimal entanglement 
witness based on second moments, with respect to a state with co- 
variance matrix 7, as a separating hyperplane from the convex set 
of second moments consistent with separable Gaussian states (dark 
grey), (b) Curved quadratic witness, which is curved towards the 
above convex set. 



Commonly used tests for entanglement are entanglement 
witnesses jj^. An entanglement witness is an Hermitian oper- 
ator W such that Tr[l/l^as] > holds for all separable states 
(Is, while Tr[FF/9] < for some entangled state p L2J. These 
tests are based on expectation values. In contrast, the tests 
mentioned above correspond to entanglement witnesses based 
on second moments, or equivalently, on 'variances' or 'uncer- 
tainties'. The covariance matrices collecting the second mo- 
ments of the states that lead to a given fixed value on the left 
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hand side of Q give rise to a hyperplane in the set of all sec- 
ond moments embodying the correlations of the state, see Fig- 
ure ^ It is a separating hyperplane: all second moments that 
correspond to separable states - as well as some correspond- 
ing to entangled states - are on one side of the hyperplane. 
Hence the test either confirms the presence of entanglement 
or returns an inconclusive result. 

These criteria are different than, for example, criteria di- 
rectly based on the positivity of the partial transpose (PPT) 
fW\, requiring the full knowlegde of all second moments. 
This PPT condition then also constrains symplectic eigenval- 
ues of covariance matrices, and hence marginal and full puri- 
ties 103 El- 

It has been shown how the set of entanglement witnesses in 
this sense can indeed be completely characterized fTT]. How- 
ever, important open questions remain: Given a state with 
some covariance matrix, can one find an optimal entangle- 
ment witness in the sense that it most robustly detects the 
state as being entangled? This is important, since for exper- 
iments aiming at the production of a specific entangled state, 
the answer would deliver an optimal test detecting the en- 
tanglement, optimally robust to noise. Another question is 
whether this can be done for all different separability classes 
versus multi-particle entanglement. This paper gives a pos- 
itive answer to these questions. In particular, we show that 
all such optimal tests (and not only the test whether a covari- 
ance matrix corresponds to a separable Gaussian state 1 14]) in 
the bi- and multi-partite setting arise from solutions to certain 
semi-definite problems fT3\ . 

To provide a service to the community, we present a pub- 
licly available software package, consisting of the functions 
FullyWit and MultiWit. Given the covariance matrix of 
a state, the first one finds optimal witnesses detecting entan- 
glement when a certain splitting of the parties sharing the state 
is held fixed, while the latter identifies witnesses detecting 
only genuine multi-partite entanglement. Further, the routines 
allow to restrict the types of measurements that one wants to 
perform. For example, for a two-mode squeezed state in the 
so-called standard basis under no further constraints the test 
as in Eq. Q would be delivered as output. Hence all entan- 
glement witnesses in the bi- and the multi-partite setup are 
efficiently obtained. 

Further, we show that for each entanglement witness, there 
is a one-to-one correspondence to a curved quadratic witness, 
thereby generalizing the known product criteria \ 16], related 
to expressions of the type 

{(u-{u)pf)p.{{v-(v)p)^)p<]. (2) 



definite programs designed to find the optimal witnesses de- 
scribed above are presented in the Sections IIVI and |V] We 
make some remarks on the form of the witnesses when the 
state in question does not include correlations between posi- 
tion and momentum variables and other issues in Section lVll 
The theoretical part is concluded with the characterization of 
the product witnesses in Section IVTIl The final Section IVlIII 
contains several numerical examples where optimal witnesses 
have been found for several states with the help of the func- 
tions FullyWit and MultiWit based on the results of the 
Sections lIVl andlVl respectively. 



11. PRELIMINARIES 

A. Definitions 

We consider system consisting of n modes, associated with 
canonical coordinates r = {xi ,Pi,X2,P2i ■ ■ ■ , , Pn), satis- 
fying the canonical commutation relations, giving rise to the 
familiar skew symmetric matrix 
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These canonical coordinates will typically - but not neces- 
sarily - refer to amplitude and phase quadratures of a finite 
number of modes of the electromagnetical field of light. In 
the following we will often consider the multi-partite situ- 
ation: here, a subsystem A embodies ha modes, system B 
consists of TiB modes, C of nc modes, and so on, such that 
TiA + ns + ■■■ — n. In this paper, we will investigate entangle- 
ment and separability properties of such states on multi-partite 
systems. We will refer to a split as a coarse graining of sub- 
systems, i.e., a distribution of the n physical subsystems into 
groups that are considered the subsystems. 

For our purposes the moments of the states will play the 
central role. The first moments are the displacements in phase 



space, dj 



0, ] 



1 ..... n. The second moments, the 



variances, can be collected in the covariance matrix 7 of the 
state, with entries 



7j,fe = 25R((fj - {fj)p)(fk - {fk)p))p 



(4) 



J, fc = 1, .^_^n. For a survey about these preliminaries, see 
also Refs. 1171 fl81. Any covariance matrix of a quantum state 
satisfies the Heisenberg uncertainty principle 



These correspond to separating hyperplanes which are curved 
towards the set of separable covariance matrices, cf. Figure[2 
Also, the complete set of generalized quadratic witnesses is 
stated and discussed. 

The paper is organized as follows: In Section|II| we recall 
the basic definitions regarding Gaussian states, the classifica- 
tion of entangled states in a general setting, and semi-definite 
programs. In SectionUll] we state the definition of linear en- 
tanglement witnesses based on second moments. The semi- 
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in turn, for every real symmetric matrix 7 e R' 
isfying (jSj there exists a physical state with just these sec- 
ond moments 7. Gaussian states - those quantum states for 
which the characteristic function is a Gaussian in phase space 
- are uniquely characterized by their first and second mo- 
ments ilTt II8II . Such Gaussian states play a central role 
in continuous-variable quantum information, essentially since 
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the Gaussian operations - completely positive maps preserv- 
ing the Gaussian character - are to a large extent readily ac- 
cessible However, the criteria we present can also 
detect the entanglement of non-Gaussian states. 

It is a very useful fact that any matrix M e ]R,2mx2m^ ^ 
0, can be diagonahzed as 

SMS'^ ^ D, (6) 

where S e Sp{2ni, R) is not an orthogonal matrix, but one 
leaving the symplectic form invariant, i.e., SaS^ = a. These 
canonical transformations are referred to as symplectic trans- 
formations. The diagonal matrix D can be taken to have the 
form D ~ (si, si, . . . , s™, Sm) with si, . . . , > 0. These 
values are the symplectic eigenvalues of M (different from 
the eigenvalues), which are also given by the eigenvalues of 
the matrix A/^/^(i(7)Af The symbol str denotes the sym- 
plectic trace of a matrix which is defined as 

m 

str[A/]=^s™, (7) 

counting each symplectic eigenvalue only once. So it is essen- 
tially the trace of the matrix after symplectic diagonalization. 

B. Separability 

Quantum states of bi-partite systems may be classically cor- 
related or entangled. If they can be prepared by means of lo- 
cal quantum operations and shared randomness alone, a state 
is called separable. A state vector 10) in turn is entangled if it 
cannot be written as as a tensor product 1 0) = \a)®\h) of state 
vectors |a), \b). A gene ral mixed state of a bi-partite system 
is called separable flY] if it can be represented as a convex 
combination of products 

^P#.>('/'.l®IV^.)(^.|, (8) 

i 

where pi > and J^iPi = 1- Otherwise, it is called entan- 
gled. 

For quantum systems consisting of more than two con- 
stituents, different kinds of classically correlated, i.e., of sep- 
arable states are conceivable. Depending on possible prepa- 
ration strategies, a state can be classified according to a full 
hierarchy. At the lowest level of the hierarchy are those states 
that contain no entanglement at all. Such an TV-partite (mixed) 
state p is called /M/fy separable, if it can be written as a convex 
combination of product states, so as 

P = E^'fclV'')^^^!^'^ ® I'A.X'^.I^'^ (9) 

k 

where Pi > for all i and 'Y^^Pi = 1. These are states that 
can be prepared by means of local operations with respect to 
all subsystems, together with shared randomness. 

In a multi-partite system, yet, also other classes of sepa- 
rability are possible. To obtain a hierarchy, one may consider 



fc-partite splits, where each of the parts is considered a subsys- 
tem in its own right, and refer to states that are fully separa- 
ble with respect to such a fc-partite split as bein^ k-separable 
ll22l (for a short review, see also Ref. [24]). Towards 
the end of this hierarchy are the 2-separable or bi-separable 
states: they are those states for which there exists a bi-partite 
split such that the state is separable with respect to this split. 
Needless to say, a state that is fully separable with respect to 
one fc-partite split might be entangled when another fc-partite 
split is considered. Hence, all possible splits for all possible k 
have to be considered for a complete classification. This clas- 
sification is treated theoretically in Section llVl and practically 
by the routine FullyWit. 

According to this classification, a state is genuinely A^- 
partite entangled if it is not separable with respect to any split. 
However, in general, there exist states which can be written as 
a convex combination of certain fc-separable states which are 
not separable with respect to any split. This motivates the al- 
ternative definition that an A^-partite state is called genuinely 
A^-partite entangled if it cannot be written as a convex com- 
bination of some fc-separable states fzs"] for any fc > 2. For 
example, in a tri-partite system consisting of parts A, B, and 
C only those states are then genuinely three-partite entangled 
which cannot be written in the form 

PBS = Ai^p(^I^^Vfe>(^fel'^'®l0fe>(0fel^'''''(lO) 
fc 

+ A2^pf^l^)|^,)(V.,|(-^^)0|<^,)(0,|(^) (11) 
j 

where X, p'-^\^'^'\ p'-^^^'^K andp^"**^!^) form probability dis- 
tributions. Genuine A^-partite entanglement according to this 
definition is treated theoretically in Section|V]and practically 
by the routine MultiWit. 

The definitions from above immediately carry over to con- 
tinuous variable systems with canonical coordinates. Here, 
they can be expressed in terms of covariance matrices. Let 
7 be the covariance matrix of a state on n modes with finite 
second moments, which is fully separable with respect to n 
subsystems. Then there exist covariance matrices 7*-*ji = 
1, . . . ,n, corresponding to the n subsystems, such that I2ql 

7>7(^)® ...©7^"^ (13) 

Conversely, if this holds, then Gaussian states with the co- 
variance matrix 7 are separable. Hence if this criterion is vi- 
olated, then the corresponing state is entangled, irrespective 
of whether it is Gaussian or not. If it is not violated, then a 
Gaussian state is separable while a non-Gaussian state might 
be entangled. Keeping this in mind we will call all covari- 
ance matrices fulfilling Eq. il3\ fully separable to keep our 
notation simple. 

Note that the problem of testing whether il3\ can be sat- 
isfied is a semi-definite problem in its own right fl4il . see 
the subsequent subsection. It is a feasibility problem, so 
the question is whether or not such matrices 7'-^'', . . . ,7''"-' 
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can be found satisfying in turn the semi-definite constraints 
7^-') + ia > for all j = 1, . . . , n. 

This statement can be generalized to bi-separable states in 
the sense of a convex combination of pure bi-separable states. 
Let 7bs be the covariance matrix of a bi-separable n-partite 
state with finite second moments. Then there exist partitions 
TT of the n modes into two subsystems consisting of to < n 
modes and n — m modes, covariance matrices 77r(fc) which 
are block diagonal with respect to the partition 7r(fc), and a 
probability distribution A so that 

7BS - ^kl^(k) > 0. (14) 

k 

Conversely, if this holds, then Gaussian states with the co- 
variance matrix 7bs are bi-separable. The same comment on 
non-Gaussian states that we made after the criterion for full 
separability applies here which we have to keep in mind when 
we call all covariance matrices fulfilling Eq. (I14> bi-separable. 
Note further that these states include those with are separable 
for a split with fc > 2, hence all states which are not in this 
set are genuinely A^-partite entangled according to the second 
definition. 



C. Semi-deflnite problems 

Semi-definite programs (SDP) [15] will play a prominent 
role in the subsequent argument. Readers familiar with the 
subject may wish to jump to the next section. Semi-definite 
programs are convex optimization problems |27] of a specific 
form: one minimizes a linear function, subject to a semi- 
definite constraint. Many problems in quantum information 
science can be cast into this form |43], essentially originating 
from the fact that semi-definite constraints appear in condi- 
tions to quantum states, as well as to quantum operations via 
the duality between positive operators and completely positive 
maps. Also, even global optimization problems can be relaxed 
to semi-definite form, with several applications to quantum in- 
formation problems ll44ll45ll . 

More specifically, a semi-definite program (SDP) is an op- 
timization problem of the following kind: 

minimizea; c^x (15) 

t 

subject to F{x) Fo + ^ FiX^ > 0, 

4=1 

where the minimization is performed with respect to a real 
vector X of length t. The problem is specified by the vector 
c € R* and the Hermitian matrices Fi £ (C^'' ,i = 1, . . . ,t. 
This is the form that is usually referred to as being the primal 
problem. It is desirable to formulate a given problem as an 
SDP, not the least because these can be solved efficiently, for 
instance by using interior point methods fT5] . 

Via Lagrange-duality the Lagrange-dual of the above prob- 
lem can be formulated. Dual problems to SDPs are again 
SDPs, where essentially the roles of objective variables and 



constraints are interchanged. The so-called dual problem can 
be formulated as follows 

maximize^ — Tr[Fo^] (16) 
subject to Z > 0, 

Tt[F,Z] = c. 

The objective value of every solution of the dual provides a 
lower bound to the value of any solution to the primal problem 
and vice versa. This is referred to as weak duality: For feasible 
X and Z, i.e. x and Z fulfilling the respective constraints, 

Jx + Tt[FoZ] = Tt[F{x)Z] > (17) 

holds, where the inequality is due to the fact that F{x) > 
and Z > 0. If either the primal or the dual problem (or both) 
are strictly feasible, meaning that there exists a feasible vector 
X such that F{x) > or there exists a feasible Z > 0, then 
there exist a;* and Z* such that 

c^x* = -Tr[i^o^1- (18) 

This is referred to as strong duality. 

An important class of problems are the feasibility problems. 
Here, c 0, so that the primal problem amounts to checking 
whether there exists any feasible x fulfilling the primal con- 
straints. In this case, Tr[FoZ] > has to hold for all feasible 
Z. Hence, if there is a feasible Z with Tr[Fo2'] < 0, then the 
primal problem cannot be feasible. 

III. LINEAR ENTANGLEMENT WITNESSES BASED ON 
SECOND MOMENTS 

A. Concept of entanglement witnesses based on second 
moments 

We now turn to entanglement witnesses based on second 
moments. They are tests for entanglement based on linear 
combinations of second moments. Each test corresponds to 
a hyperplane in the set of second moments. In turn, these 
hyperplanes encode the physical set-up, the type of measure- 
ment that is being physically performed. Such hyperplanes 
are defined by a real symmetric positive semi-definite (PSD) 
matrix Z and a number c e R, via the Hilbert Schmidt scalar- 
product. The hyperplane consists of all 7 such that 

Tr[Z7] = c. (19) 

So an entanglement witnesses based on second moments can 
be characterized by a real matrix Z > satisfying 1 1 1 ] 

(i) Tr[Z7,.] > 1 for all (fully) separable 7s, (20) 

(ii) Tr[Z7] < 1 for some entangled 7. (21) 

As mentioned in the introduction, this notion is very similar 
to the one for entanglement witnesses based on expectation 
values: note that here, however, the witness refers to second 
moments of quantum states. If condition (ii) is fulfilled, then 
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the quantum state giving rise to 7 is entangled, irrespective of 
whether it is Gaussian or not, as discussed in Sec. Ill Bl 

These hyperplanes are hence separating entangled covari- 
ance matrices from the set of fully separable covariance ma- 
trices. That this characterization is possible is due to the fact 
that the set of fully separable covariance matrices is convex 
and closed. Its boundary is given by the matrices (B^^iJk 
fulfilling the condition ®^!^]^7/c > icr. It is clearly convex, 
because the semi-definite constraint defining fully separable 
covariance matrices is preserved under convex combination: 
if ffife=i7fe > *f ©fc^i% > ia then 



)Li 7fe + (1 - a) ©fe=i Vk > 10 



(22) 



for a e [0, 1]. Further, the set is closed. Firstly, the subset of 
covariance matrices of the form ®^'^]^7fc is closed itself since 
its complement is open: if a matrix 7 has nonvanishing off- 
diagonal elements, then in its neighborhood there will be only 
matrices with nonvanishing off-diagonal elements. Secondly, 
the constraint 



(23) 



defines a closed convex cone which is a subset of the space of 
matrices ©"^^-{^k- Because the set is convex and closed, there 
exist hyperplanes separating a covariance matrix 7 which is 
not separable from the set of separable covariance matrices 



B. Familiar examples 

Any such matrix Z encodes the measurement pattern per- 
formed for a certain test, so the linear combination of second 
moments that is required to assert that a state was entangled. 
For example, the matrix Z for the familiar test of Eq. ([ij from 
Ref . 01 can be written as 
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10 10 
10-1 
10 10 
0-101 



(24) 



Using the definition (0} it follows that 

Tr[Z7] = {{u - {u)pf), + {{v - (25) 

such that inequality Q is equivalent with 

Tr[Z7] < 1, (26) 

from which one can conclude that the state must have been 
entangled. In turn, the test gives rise to a hyperplane in the 
space of second moments separating a subset of all entangled 
states from the separable states. The general criterion of Ref. 
10], dependent on a parameter a e ]R\{0}, reads as 



{(Ua - {Ua)pT)p + {{Va " {Va) pf) p > + 



1 



for all separable p, where 

Ua = \a\xi + -X2, Pa^\a\pi~-P2- (28) 
a a 

This test corresponds to a particular hyperplane for each a, 
corresponding to the entanglement witness 



(29) 
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so that 



Tr[Za7s] > 1 for all separable 7s, 
Tr[Za7] < 1 for some entangled 7. 



(30) 
(31) 



(27) 



In Ref. fTl fToll it was further shown that the covariance matrix 
can always be brought into the standard form of a direct sum 
in position and momentum variables, using local symplectic 
transformations. The test ( I28t is necessary and sufficient for 
covariance matrices in this standard form. In other words, if 
the covariance matrix of a state is known, and we allow for 
further local operations based on that knowledge, then the set 
of witnesses Za completely characterizes the set of separable 
Gaussian states. 

Another example of criteria that correspond to witnesses 
based on second moment are the tests for multi-partite entan- 
glement of Ref. 1 8]. In turn, the set of witnesses of second 
moments can be fully characterized fllll . as sketched in the 
Appendix. 



C. Detecting entanglement in non-Gaussian states 

There is finally a comment in order concerning non- 
Gaussian states: As mentioned before, all such tests also de- 
tect non-Gaussian states as being entangled. This is due to 
the fact that if a Gaussian state is entangled, then every non- 
Gaussian state with the same second moments is necessarily 
entangled. Hence, one does not have to assume a priori that 
the state in question is exactly Gaussian. This knowledge is 
typically not available anyway, without referring to full quan- 
tum state tomography. 

One has to be aware, however, that in such infinite- 
dimensional quantum systems, the entangled states are trace- 
norm dense in state space |i28.3. That is, in any neighborhood 
of a separable state an entangled state can be found, and the 
volume of the set of separable states has measure zero. This 
situation remains unaltered if one introduces a constraint to 
the mean energy of the system fSs']. 

Nevertheless, it still makes sense to test for entanglement 
in this setting: one only has to state the result in the form that 
if a state is detected as being entangled then a trace-norm ball 
centered at it is not consistent with any separable state. In 
other words, one can still meaningfully say that a state is 'far 
from being separable'. 
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IV. TASK OF FINDING WITNESSES AS SEMI-DEFINITE 
PROBLEM 

The common situation that one encounters is the following: 
one knows what kind of entanglement one would like to see 
in a certain state, prepared in some setup. Also, one typically 
has an idea about how the covariance matrix 7 of the prepared 
state roughly looks like or at least about how it is desired to 
look like. The question is: what measurements have to be 
performed in order to most easily detect the entanglement? 
More specifically, in case of bi-partite entanglement, the task 
is the following: 

We have two parts consisting of ua and n b modes, respec- 
tively. For a given 7 (the covariance matrix that we suspect 
that we have) corresponding to an entangled Gaussian state 
we would like to find the test Z with the property that Z de- 
tects 7 as corresponding to an entangled state, such that 

w = Tr[Z7] < 1 (32) 

takes its minimal value. This is the test which 'most distinctly' 
detects 7 as originating from an entangled state, in a way that 
is most robust against detection imperfections. Geometrically, 
we aim at finding the hyperplane with the greatest distance 
from 7. Obviously, not only 7 is detected as coming from an 
entangled state by this test, but the test is optimized for this 
specific guess. 



A. The primal problem 

It turns out that the previous problem is related to the fol- 
lowing optimization problem. For a separable 7, we have 
that there exist covariance matrices ja and 7s, satisfying the 
Heisenberguncertainty relation 7^ 07b + i(7 > 0. For covari- 
ance matrices 7 corresponding to entangled states, we may 
write the primal problem in the following form: 

minimize^^^^^,^^ {-Xe), (33) 
subject to 7 — 7A ffi 7b > 0, 

©7i3 + il+Xe)i(T > 0. 

If there is an optimal solution with Xe > 0, then 7 is separable, 
because 7a©7b fulfils an even stricter form of the Heisenberg 
uncertainty relations. If Xe < 0, then 7 is entangled, since 
7/1 ffi 7_B can now violate the uncertainty relations. This is 
actually just the p-measure from Refs. 12011 . uptop = 1/(1 + 
Xe). For Gaussian states, — log2(p) is a lower bound for the 
logarithmic negativity, defined as 

En{p)= log^Wp^Wi, (34) 

where denotes the partial transpose of p and ||.|ji is the 
trace norm. It is moreover identical to the logarithmic nega- 
tivity for 1 X n-mode systems 1I20I1 . The negativity ll29il is a 
measure of entanglement, and indeed a monotone under local 
operations and classical communication L30..3L»32|] . 



B. The dual problem 

The dual problem of the above problem can easily be found. 
The key point in what follows is that from the dual of the 
above semi-definite problem allowing for a pre-factor in the 
Heisenberg uncertainty one can extract the required optimal 
tests. The dual problem can be cast into the form 

maximize^ -Tr[(7 © ia)X] (35) 
subject to X > 0, 

Tr[(0®iCT)X] -1, 
Tr[(-^^j-fc©F,-fe)X] =0, 

j,k=l,.. .,nA, 

j, fc = + 1, . . . ,n, 

where the maximization is performed over Hermitian matrices 
X £ (|2;4nx4n jj^j^ matrix corresponds to a partitioning of 
degrees of freedom labeled 1, . . . , of system A first, then 
ua + I, ■ ■ ■ ,n labeling the modes of system B, and then again 
the same ordering. Fj^h, j,k = 1, . . . ,n, form a set of real 
symmetric matrices all entries of which are zero, except 

iF,,k)j,k = {F,.k)k^, - 1. (36) 

These matrices Fj^k form a basis of all real symmetric n x n- 
matrices. The form of Eq. (I35> of the dual problem becomes 
manifest when expressing the primal problem in terms of this 
operator basis, and writing the two semi-definite constraints of 
the primal problem in form of a direct sum as one constraint. 

Due to the block diagonal structure of (7 © ia) and all con- 
straints, we can without loss of generality assume that 

X^Xi®X2. (37) 

The first constraint is then equivalent to Xi 2 > 0. The 
latter constraints in the dual problem lead to Tr[Fj.fcXi] = 
Ti[Fj^kX2], which restricts the real symmetric single system 
blocks of system A and B in Xi and X2 to be equal. This 
generalizes directly to the case of N subsystems: then the N 
real symmetric single system blocks of Xi have to be equal to 
the N real symmetric single party blocks of X2- The matrix 
X2 is further restricted by the condition 

Tr[(0 © ia)X] = Tr[icrX2] = -1. (38) 

Finally, the dual objective function is 

Tr[(7©icr)X] ^ Tr[-fXi] + Ti[iaX2] = Tr[7Xi]-l. (39) 

Since 7 is real and symmetric, we have further that Tr[7Xi] = 
Tr[7Xf ], where Xf is the real part of Xi. 

To summarize, the dual problem can be formulated as 

minimizexi.x^ T^iilXf] - I, (40) 
subject to = X^'^'"^ 

Xi >0, X2> 0, 
Tr[iaX2] = -1. 
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In this formulation, no basis is used explicitly. XJ' 2 refers 
to the block diagonal matrices obtained from through 
pinching to the blocks of system A and B (i.e., a projection 
onto the form of a direct sum). Now we are in the position 
to formulate the connection between the separability problem 
and witnesses based on second moments: 

Proposition 1 (Optimal witnesses in bi-partite systems) 

For every feasible solution X to the dual program formulated 
above, the matrix X^" fulfils the witness condition i20i . If 7 
is entangled, then 

Tr^Xf] < 1, (41) 

so that X I also fulfils condition i21l . Further, Tr[7X[''] is the 
minimal value of Tr^y Z] for any witness Z. 

Proof. From weak Lagrange duality ( I17> it follows that 

c^x + Ti[{j(Sia)X]>0, (42) 

where c is the vector specifying the objective function of the 
primal problem, being equivalent with 

TT[-fX'^]>l+Xe. (43) 

In this case, there is always a strictly feasible X: just take 
X = 1 © (1 — ia/{2nA + ins))- Hence there exist feasi- 
ble X and X such that equality is obtained in Eq. ( I43> . We 
have seen before that a:e > for all separable states. Hence 
the condition ( I20> is fulfilled. On the other hand, if 7 is an 
entangled covariance matrix, then Xe < 0, hence also the con- 
dition (12 1> if fulfilled. Since the equaUty holds in Eq. j43> , 
Tr[7X]^''] reaches the minimal value. q 

An analogous proposition holds for witnesses detecting full 
separability in a system of N subsystems. Hence all classes 
of /c-separability can be tested with this criterion. 

V. DETECTING GENUINE MULTI-PARTITE 
ENTANGLED STATES 

The previous section was devoted to tests of full separa- 
bility. Here, we formulate the problem for excluding bi- 
separability for systems of N subsystems, each consisting of 
Ui modes, i = 1, iV. As before, the total number of modes 
is n — J2i ^i- The condition bi-separable states have to fulfil 
is inequality ilAl . 

A. Primal problem 

We write the primal problem in the following form: 

minimize{^^(^j}^^^ -Xe, (44) 
subject to 7 - X! ^^(■'^I - ^ 

k 

lTr{k) + ^feicr > for all A:, 

Afc = 1 + Xe, 

fe 

Afc > for all k. 



Here 7r(fc) are all K = 2^~^ — 1 possible bi-partite partitions 
of the N systems. The matrices 7^(fe) are block diagonal with 
respect to the partition 7r(fc). If the solution Xe > 0, then 7 
is bi-separable, because the matrices 77r(fc) fulfil an even 
stricter form of the Heisenberg uncertainty relations. If the 
solution Xe < 0, then 7 is genuinely multi-partite entangled, 
since 77r(fc) can now violate the uncertainty relations. 

With the basis introduced in the previous section, the prob- 
lem can be formulated as 

minimize^^x(fc)j -Xe, (45) 

bd,re,7r(fc) 

subject to 7+ ^ (-i^ijXf^>0, 

bd,re,7r(/c) 

^ ^ + Afcicr > for all k, 

Afc - Xe - 1 > 0, 

fc 

-(I]Afc)+Xe + l>0, 

Afc > for all fc, 

where the index 'bd,re,7r(fc)' refers to 'block-diagonal and real 
with respect to the partition 7r(fc)'. The set of constraints can 
again be cast into the form of a single constraint in a direct 
sum form. This helps to identify the respective terms in the 
dual problem. 

B. Dual problem 

We can again assume the Hermitian matrix X E 

(^(2n[K+i]+2+K)xi2n[K+i]+2+K) to be block diagonal. The 

dual problem is then given by 

maximize^ -Tr[(7 02„k ® (-I1) li 0^^)^], (46) 
subject to = x'^MMk) ^jj 

TT[iaXk+i] + Xk+2 - Xk+3 + Xx+z+k = 
for all fc, 

Xk+2 — Xk+3 — 1, 

where l^; and 0^ are the d-dimensional identity operator and 
operator, respectively. With the last constraint, the objective 
function reduces to 

Tr[(7 ® Q2nK ® (-I1) ® li ® Qk)X] 
= Tr[7Xr] - {Xk+2 - Xk+s) 
= Ti[-iX'^] - 1. (47) 

Now we can formulate the connection between the bi- 
separability problem and witnesses: 

Proposition 2 (Witnesses for multi-partite entanglement) 

For every feasible solution X to the dual program formulated 
above, the matrix Xi satisfies the witness condition 

Tt[X['-ibs] > 1 (48) 
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for all bi-sepamble "/bs- If ^ is genuinely multi-partite entan- 
gled, then Tr[7X['^] < 1 so that X" also satisfies condition 
MIX . Further, TrYiX'i] is the minimal value of Ti^yZ] for 
any witness Z detecting only genuinely multi-party entangled 
states. 

Proof. From weak Lagrange duality Mil it follows that 

c^a; + Tr[(7e02„x©(-li)©li©0K)X] > 0,(49) 
which is equivalent with 

Tr[7Xr] > l + xe. (50) 
In this case, there is always a strictly feasible X: just take 

1 



Xk+2 



1, X, 



fe+1 



1 



2' 

XK+3+k — 1 



X 



K+3 



(51) 

(52) 
(53) 



for all k. Here, n is the total number of modes as before. 
Hence there exist feasible X and x such that equality is ob- 
tained in the last equation. We have seen before that > 
for bi-separable states. Hence fulfils the condition ilOi . 
Further, if 7 is genuinely multi-partite entangled, then x^ < 0. 
Hence also respects the condition i2l\ . Finally, since 
the equality holds in Eq. i5Q\ . TT[-fXf] reaches the minimal 

□ 



value. 



VI. REMARKS 

A. Direct sum form of tests for block diagonal covariance 
matrices 

Often, the covariance matrices of generated states exhibit a 
direct sum form with respect to position and momentum vari- 
ables. In the simplest two-mode form, this corresponds to a 
covariance matrix of the form 



Ci ^5 

6 

C5 6 

^6 C4 



with some ^1 , . . . , € Not only that every two-mode co- 
variance matrix can be brought into such a form by means of 
appropriate local symplectic transformations, but many of the 
second moments of typical generated states exhibit approxi- 
mately this form anyway. Two-mode squeezed states or noisy 
variants have covariance matrices of this form, needless to say. 
Then, the question is: can the test also be taken to be of this 
form, without losing optimality? The well-known test in ([0 
from Ref. to give an example, is of such a form. This 
question can be positively answered. 

Proposition 3 (Direct sum in position and momentum) 

Let ^ be a covariance matrix of the form of a direct sum 



of matrices corresponding to position and momentum coor- 
dinates (as in Eq. i l54t for two modes). If Z is the optimal 
solution of the dual problem for 7, then the pinching 



where 



Z' — PpZPp + PxZPx, 



Pp = diag{0, 1, ... ,0,1), 
Px = diag{l,0,...,l,0) 



(55) 

(56) 
(57) 



is also a feasible solution with the same objective value 
Tr[jZ]. 

Proof. The first steps of the proof are straightforward: Since 7 
is of the form of a direct sum of matrices corresponding to the 
position and momentum coordinates, 7 ~ PxjPx + PplPp 
holds, and hence Tr[Z7] = Tr[Z'7]. Now it has to be shown 
that Z' is a witness. First, it fulfils Z' > 0, as every prin- 
cipal submatrix of a positive matrix is positive. Secondly, 
Tr[Z'7s] > 1 has to hold for all separable 7^. This is equiva- 
lent to Tr[Z'7s] = Tr[Z7^] > 1, where 



7.- 



PxlsPx + PplsPp 



(58) 



Hence if 7^ is a separable covariance matrix then Z' is a wit- 
ness. The covariance matrix 7s fulfils 7s — 7yi ffi 7b > for 
some covariance matrices 7^,75. Then, clearly 7s — 7^ 
7b > holds. It remains to show that 7^ ® 7^ + ^o" > 
if 7A © 7s + icr > 0. Due to the block-diagonal struc- 
ture, it suffices to show that for any covariance ?/ satisfying 
ry > id, also 77' > ia holds. It is in this context conve- 
nient to order coordinates as {xi, . . . , af„,pi, . . . ,Pn)- Then, 
T]' is obtained from rj as a result of a pinching. In the fol- 
lowing, we use standard notation from Ref. |17]. The idea 
of the proof is that we use appropriate symplectic transfor- 
mations that commute with both P^ and Pp to transform 
PxriPx+PpVPp+if^ PxMPx + PpMPp+ia, such that the 
problem is reduced to a single mode problem. First, the posi- 
tion part of 77 can be brought to diagonal form by the congru- 
ence 7 (O © 0)7(0 ® O)^, where O G 0(n). Therefore, 
[O, Px] — [O, Pp] — 0. Then, with single mode squeezings. 



S = diag(di, ...,dn, 1/di, . . . , l/rf„) 



(59) 



di e IR\{0} for 1 = 1, ...,71, the position part can be made 



^^^^ proportional to the identity. Again, [S, Px] — [S, Pj 



0. 

Finally, the momentum part can be made diagonal, using an 
appropriate V (S) V, V E 0{n), again of the form of O from 
the first step, leaving the upper block invariant. Hence, 

M =iV® V)S{0 ® 0)?7(0 © OfS^iV © Vf (60) 

is diagonal in both the position and the momentum part. 

We can apply a pinching such that also the off-diagonal part 
of M is diagonal, leaving a invariant. Since M + ia > 
holds, then M' + ia > holds as well, where M' is the 
pinched form of M. The covariance matrix is now a direct 
sum of single modes. But then, the validity of the statement 
becomes obvious: if 



a c 
c b 



+ ia>0. 



(61) 
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then always also 



a 
b 



ia> 



holds true. Hence, finally we arrive at 7^ ® 7^ > icr. 



(62) 



□ 



Hence, it does not restrict generality for covariance matri- 
ces 7 in form of a direct sum of position and momentum con- 
tributions to take a test with the same form. In any case, even 
if 7 does not have this form, one may look at such tests, which 
we will later see again in the context of product criteria. 



B. Incorporating practical measurement constraints 

Often, some combinations of second moments are more ac- 
cessible via experiment than others. One is typically well- 
advised to avoid estimating all entries of the covariance ma- 
trix, but only directly those combinations that are required 
- as is, e.g., routinely done using the above test Q. In the 
program, to be described later, additional constraints to incor- 
porate specifically accessible measurement types (for exam- 
ple via interferometers, rather than through homodyning mea- 
surements) can be taken into account via a finite number of 
linear constraints of the form 



Tr[i?,Z] = 0, 



(63) 



i ~ 1 , . . . , /. If this set of / constraints is too restrictive, it 
may happen that no test is found that can detect the entangle- 
ment. 



C. Remarks on quantitative statements 

One should be tempted to think that whenever a state 'vio- 
lates' such a criterion by a large degree, it should be very much 
entangled, in quantitative terms. For Gaussian states this is a 
simple issue, as the result of the test is essentially just a lower 
bound to the logarithmic negativity, see above (compare also 
Ref. |33]). More relevant, yet, are statements that do not as- 
sume the Gaussian character of the state, as the very point of 
the test is that one not only does not need full tomographic 
knowledge of the quantum state, but not even knowledge of 
all of its second moments. 

Yet, it is not true that the Gaussian state is the one with 
the smallest logarithmic negativity, given some value in Eq. 
Q: there can be small violations. It has been shown that 
for a given value in Eq. Q, non-Gaussian states may have a 
slightly smaller logarithmic negativity [34]. Another exam- 
ple has been presented in Ref. I3al . where it has been shown 
that for a fixed full covariance matrix (and not only a fixed 
value of Eq. Q), there exists a non-Gaussian state assuming 
a smaller negativity than the corresponding Gaussian state. It 
is still true, however, that for a given left hand side of Eq. Q, 
one can still find a lower bound of the logarithmic negativity, 
indicating that 'a state that very much violates this criterion is 
also very much entangled' . 



For the entanglement of formation, the smallest degree of 
entanglement is in turn assumed for a Gaussian state for sym- 
metric 1 X 1-mode states |35], as well as for the squashed 
entanglement l36ll . Also, the conditional entropy. 



C{p) - S{pa) - S{p) 



(64) 



for states p, a lower bound to the distillable entanglement, 
takes in general its smallest value for Gaussian states i37ll . 
Hence if the full covariance matrix of a state is known, then it 
is possible to evalute the entropies of the corresponding Gaus- 
sian state with that covariance matrix, yielding a lower bound 
to the distillable entanglement of the non-Gaussian state. One 
hence does not have to assume the Gaussian character then: 
One can then measures the moments, evaluates the entropies 
of the respective Gaussian state, and take this value as lower 
bound to the distillable entanglement of the true non-Gaussian 
state. 



VII. CURVED WITNESSES: ALL PRODUCT CRITERIA 

In this section, we now turn to curved witnesses, tests that 
do not correspond to linear combinations of second moments, 
but to quadratic ones. Interestingly, the measurements that 
have to be performed are just the same ones as in linear tests, 
only the combination of the respective outcomes is different. 
It turns out that the use of quadratic tests is always advanta- 
geous to linear tests. Frankly, it is always 'better to multiply 
the numbers instead of adding them', see Figure 1. 

Geometrically, such tests - generalized product criteria - 
correspond to curved surfaces, which are not hyperplanes. 
They are curved towards the set of second moments of sep- 
arable Gaussian states. In the setting before, any witness Z 
can be decomposed as 



Z — Zx Zp 



where 



PpZpPp 

Px Px 
Pp Pp 



Px ^x.pPx 
Pp^x^pPp 



0, 
0. 



(65) 



(66) 
(67) 
(68) 
(69) 



The projectors and Pp are defined in Eq. i56\ and i51l . 
respectively. Let us first look at the class of witnesses Z that 
are of the form 



Zn 



(70) 



From each witness of this form, a product criterion of the type 
of Eq. (|2ji can be derived, which is stronger than the witness. 
Moreover, there is a one-to-one correspondence between the 
witness and the product criterion. Hence, all product tests in 
the strict sense, involving the product of variances with re- 
spect to position and momentum coordinates, are obtained. 
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Proposition 4 (All product criteria) IfZ^-^Zp is a witness, 
then 

(i) Tr[Za;7s] Tr[Zp7s] > i for all separable ^s,(Ti) 

(ii) Tr[Z2;7] Tr[Zp7] < — for some entangled j, (72) 

and the set of entangled covariance matrices detected by this 
quadratic test is strictly larger than that detected by the lin- 
ear witness + Z^. In turn, if Z^ and Zp are symmetric 
matrices fulfilling and 
the conditions (i) and (ii), then Z^ + Zp is an entanglement 
witness. 

Proof. This proof makes use of the material presented in the 
Appendix. If the witness Z^ + Zp detects a covariance matrix 
7, Th:[{Z^ + Zp)j] < 1, then clearly also Tr[Z^7] Tr[Zp-/] < 
1/4 holds. However, if Tt[{Z.j: + Zp)-/] > 1, then it does 
not follow directly that Tt[Z^j] Tr[Zp7] > 1/4. Therefore, 
we continue by showing that if Tr[Zxj] Tr[Zp7] < 1/4 then 
there exists a witness Z' such that Tr[Z'7] < 1. It follows that 
Tr[Zj^7] Tr[Zp7] < 1/4 for entangled covariance matrices 7 
only, concluding the proof of the statements (i) and (ii). 

We define Z'^ = aZ^ and Z'p = Z^/a, a e R\{0}, such 
that 

TrK7] -Tr[Z;7]. (73) 

This, together with Tr[Z^7] Tr[Z^7] = Tr[Z^7] Tr[Zp7] < 
1/4 implies that Tr[(Z^ + Zp)j] < 1. Further, the matrix 
Z' = Z'^ + Z'p fulfils the witness condition J2k ^^^'k > 1/2 
from the Appendix, where Z'^ is the block on the diagonal of 
Z' of party k: since Z'^ = SZkS"^, where S £ Sp{2nk, R) is 
the symplectic transformation 

S ^ iV^P^+Pp/V^), (74) 

the symplectic trace of Zj^ is equal to that of Zk- Hence Z' 
is a proper witness detecting 7. Note that this argument is 
independent of whether Z is a witness against full separability 
or for multi-party entanglement. 

Finally, the product criterion is detecting strictly more 
entangled covariance matrices than the witness Z^ + Zp 
since only the implication Tr[(Z2; + Zp)'j] < 1 
Tr[Z2:7] Tr[Zp7] < 1/4 holds, while the converse direction 
does not necessarily hold for some entangled 7. For instance, 
consider a witness Z^ + Zp and a covariance matrix 7 such 
that 

Tr[(Z, + Zp)7] < 1. (75) 

The covariance matrices S-fS"^, where S = {y/aPx+Pp/ y/a) 
as above, are proper covariance matrices for all a > 0. How- 
ever, not all of them are detected by Z^ + Zp, since it is always 
possible to choose an a > such that Ti[{Zx + Zp)S^S^] = 
aTr[Zj,7] + Tr[Zp7]/a > 1. In contrast, the product witness 
detects the whole family since 

Tr[ZxSjS^] Tr[Zp57^^] = Tr[Z,7] Tr[Zp7] < 1/4. 

(76) 



Therefore, the product witnesses are stronger tests than the 
respective linear tests. 

Finally, assume that there exist symmetric matrices Z^ — 
PxZxPx and Zp = PpZpPp, such that the conditions (i) and 
(ii) are fulfilled. From condition (i) it follows directly that 
tr[{Zx + Zp)js] > 1 for separable 7^. Then choose a 7 
such that (ii) holds. As before, TT[Zxj{a)] Tr[Zp7(a)] = 
Tr[Z:,7] Tr[Zp7] < 1/4 for 7(a) = SjS^, where S is the 
symplectic transformation J74> . If we pick a £ R\{0} such 
that 

Tr[Z,7(a)] = Tr[Zp7(a)], (77) 

then it follows from Tr[Z^7(a)] Tr[Zp7(a)] < 1/4 that 
Tr[{Zx + Zp)'y{a)] < 1. Hence Z^ + Zp is an entanglement 
witness. q 

As a matter of fact, if we allow for Zx.p 7^ 0, then we 
get a one-to-one correspondence of witnesses and generalized 
product criteria, fully characterizing the convex set of separa- 
ble second moments. 

Proposition 5 (All generalized product criteria) If Z = 

Zx + Zp + Zx^p is a witness, then 

(0 Pz{js) ^ ^ for all separable 'js, (78) 
( ii) Pz (7) < — for some entangled 7, (79) 

where 

Pzh) - Tr[Z,7] Tr[Zp7] 

+ iTr[Z,,p7] - i(Tr[Z,^p7])2, (80) 

and strictly more entangled covariances are detected than by 
the witness Z. In turn, if Zx, Zp, and Zx^p are symmetric 
matrices fulfilling Zx — PxZxPx, Zp — PpZpPp, and Zx^p — 
Px^x.pPp + PpZx.pPx o-nd the conditions (i) and (ii), then 
Zx + Zp + Zx^p is an entanglement witness. 

Proof. If the witness Z detects a covariance matrix 7, 

Tt[{Zx + Zp + Zx.p)'j] < 1, then also 

Tr[Z,7] Tr[Zp7] < AV4 (81) 

holds, where 

A= l-Tr[Z,,p7], (82) 

which is equivalent to Pz (7) < 1/ 4. All the other steps of 
the proof of proposition^can be performed in analogy, using 
that SZS^ = Z'^ + Z^ + Zx^p, where S = y/EPx + Pp/^/E 
as above. q 

These criteria hence form a complete set of criteria, and all 
what has been said before is also applicable to these tests. In 
a sense, these curved tests compensate for local squeezings, 
operations under which the linear tests are not invariant. Note 
also that the tests are quadratic, but still in entries of the canon- 
ical coordinates. One can also think of tests where the observ- 
ables themselves include higher polynomials. First interesting 
steps in this direction have been undertaken, for example, in 
Refs. |38i. 
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VIII. NUMERICAL EXAMPLES 



covariance matrix 7ww, detected by the witness 



We implemented the dual programs for witnesses detecting 
entanglement and genuine multi-partite entanglement in Mat- 
lab (Version 7). The routines have been made freely available 
H^]. They make use of the solver SeDuMi |40] and the inter- 
face Yalmip filll . which are also freely available. 



A. Testing full separability numerically 

The function FullyWit implements the dual program of 
Eq. ( I40t . It is called by the line 

[c Z ] =FullyWit (gamma, n, constraints ) . (83) 

The inputs are the covariance matrix gamma and a vector n, 
which holds the number of modes that each of the parties have. 
For instance, if gamma is a 6-mode state held by three par- 
ties A, B, and C, where party A holds 3 modes, party B 
holds 1 mode, and party C holds 2 modes, then n= [ 3 1 
2 ] . The symmetric covariance matrix gamma would have 
to have dimension 2n x 2n, where n = 6. Using the pa- 
rameter constraints, the witnesses can be further re- 
stricted, as explained in Section IVlll CI Until then, we will 
set constraints = 0, thereby not using this option. 

The output Z is a real symmetric matrix fulfilling the first 
witness condition Tr[Z7s] > 1 for all separable covariances 
7s. The second output is 



c = Tr[z gamma] — 1. 



(84) 



Hence if c< 0, then gamma is entangled, and Z is an optimal 
entanglement witnesses in the sense of Proposition 1 . Other- 
wise, gamma is separable. 

All the following examples were solved in fractions of a 
second on a Pentium 4 machine with 2.8 GHz and 512 Mb 
RAM by FullyWit and MultiWit, respectively. 

The first example we would like to consider is the PPT en- 
tangled state of 2 X 2 modes given in Ref. i26ll 



7ww — 
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(86) 
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(87) 



where x = 0.1394, y = 0.0374, and z ^ 0.1021. 

As a second example we consider the family of GHZ like 
Gaussian states introduced in Ref. |42]. For N modes, the 
covariance matrix is given by 

A C C ... C 
C A C ... C 



7GHZ = 



C C 



where 



and 



A 



a 
b 



C 



A C 
C A 



c 
d 



(88) 



N 
1 

N 
1 

N 
1 

N 



-2ri 



+2ri 



N - 1 

-IT' 

N ~1 



-2r2 



.+2r2 



(89) 

(90) 
(91) 
(92) 
(93) 



Further, ri > is the squeezing parameter of the first initial 
mode at the beginning of the state construction process, and 
r2 > is the squeezing parameter of the other — 1 modes 
at that stage |42]. These states have the following properties: 
they are pure, invariant under exchange of any two parties, and 
are genuinely multi-partite entangled. 

For = 3 and ri = r2 = (ln2)/2, the covariance matrix 
takes the simple form 



73GHZ = 7: 





3 


-1 


-1 





-1 


3 


-1 





-1 



-1 


3 



(94) 



The routine FullyWit yields c=-0 .500, clearly demon- 
strating the entanglement of the state for these parameters, and 
the witness is given by 



yields c=-0.1034, exemplifying the entanglement of the 



2a; -a; -a; 

2a; 2a; 2a; 

-a; 2a; -a; 

2a; 2a; 2a; 

-a; -a; 2a; 

2a; 2a; 2x 



(95) 
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where x = 0.0833 « 1/12. 

Note that the function FullyWit can be used to perform 
all the tests necessary for the first classification of multiparty 
entanglement introduced in Section HTSI For instance, sepa- 
rability across the split AB\C for three parties as in the last 
example can be tested by choosing n= [ 2 1 ] . However, if 
the split AC\B is chosen, then the parties B and C have to be 
exchanged first by transforming 7 1-^ S'^S'^, where 



10 
1 
10 



(96) 



Then the test can be performed by the command 

[c Z] =FullYWit (S'7S'^, [2 1],0). (97) 

In a similar manner, all possible partitions for N parties can 
be tested. 

B. Testing bi-separability numerically 

In analogy, the function MultiWit implements the dual 
program of Eq. j46l l. It is called by the line 

[c Z ] =MultiWit (gamma, n, constraints ) . (98) 

The inputs are again the covariance matrix gamma and a vec- 
tor n, which holds the number of modes that each of the par- 
ties have as in FullyWit. Again, we put constraint s = 
and refer to Section lvTircl 

The output Z is a real symmetric matrix fulfilling the first 
witness condition Tr[Z7s] > 1 for all bi-separable covari- 
ances7s. The second output is c= Tr[z gamma]-l. Hence if 
c< 0, then gamma is genuinely multi-partite entangled, and 
Z is an optimal entanglement witnesses in the sense of propo- 
sition 2. Otherwise, gamma is bi-separable. 

Applying this routine to the 3 mode GHZ covariance of 
Eq. ( I94> . we obtain c=-0 .305 6, showing that the state is 
genuinely multi-partite entangled. The witness detecting the 
state has the form 
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-z 
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(99) 



where x = 0.2315, y = 0.2021, z = 0.1157, and w = 
0.1875. 

Another example we consider is a four mode state occuring 
in an intermediate step of an continuous variable entangle- 
ment swapping experiment as described in Refs. i4al . Key 
steps towards full implementations and first experimental im- 
plementations have been reported recently 1 3, 47]. 

The idea can be briefly described as follows: two entan- 
gled states are produced between the modes labeled 1 and 2 



and between the modes 3 and 4, respectively. Then a cer- 
tain homodyne-type measurement is performed between the 
modes 2 and 3. In a successful implementation, after this, 
entanglement can be confirmed between the modes 1 and 4 
which have never interacted before. Here, we consider the 
state before the actual measurement, which is the state where 
the modes 2 and 3 have been brought to overlap at a 50:50 
beam splitter 

In order to see whether this can indeed be confirmed for 
realistic parameters, we construct two mode entangled states 
with a fixed degree of squeezing and a fixed degree of mixed- 
ness by 'backwards reasoning'. We start with the diagonal 
matrix 7' — diag(e^^'', e^^'', ae^' , ae^''), where r, a > 0. 
This is the desired covariance matrix after partial transposi- 
tion and symplectic diagonalization. The entanglement is re- 
flected in the symplectic eigenvalue 1 17] e^^'' < 1 for r > 0. 
We arrive at the covariance matrix 7 by the reverse transfor- 
mation, in this case we first apply a 50:50 beam splitter and 
then partial transposition with respect to the second system. 
The symplectic eigenvalues of the resulting matrix 7 are both 
equal to ^Ja. The von-Neumann entropy can for a Gaussian 
state be calculated as 1 17] 

n 

H{p) = J2 [{Nk + 1) HNk + 1) - Nk InNk] , (100) 
fc=i 

where A'^fc = (sfe — l)/2 is obtainedfrom the symplectic eigen- 
values Sfc of 7. 

For both input entangled pairs of modes, we choose r — 
2 ln(2)/3, corresponding to 4dB squeezing of the initial state 
(the degree of squeezing can be read off the smallest eigen- 
value of the covariance matrix) and a = 5, leading to H{p) — 
2.152. The covariance matrix of the state after the beam split- 
ter between the modes 2 and 3 is then given by 
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(101) 



where x = 6.4980 and y — —4.3142. The routine 
MultiWit returned c=-0.2305, demonstrating the gen- 
uine four-partite entanglement of the state. The corresponding 
witness is of the form of Eq. dlOH . where now x — 0.2352 
and y = 0.1660. 

The function FullyWit returns c=-0 . 6031 for n= [1 
1 1 1 ] . The corresponding witness is of the same form as 
Eq. ([Toll , where now X = 0.125 and y — 0.0884. However, 
this specific instance of a test can only clarify that the state is 
not fully separable. 



C. Further experimental constraints 

If the experimental set up limits the possible set of tomo- 
graphic measurements, then these constraints can be taken 
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into account by requiring that the witness Z fulfils 

Tr[ZA] = 0, 



(102) 



where A is an operator describing the measurement. We can, 
for instance, set the element Zj ^ to zero with the matrix A — 
Fj^k defined in Eq. (I36> . In this way it is possible to restrict 
Z such that no position-momentum correlation measurements 
are required, as discussed in Sec. IVl Al 

If we require that Eq. ( I102> holds, then the only constraint 
that is affected in the primal program is 



7 > 7A © 7B J + XaA> JA® JB, 



(103) 



and in analogy for the program detecting only genuine multi- 
partite entanglement. Here, xa is a new real variable in addi- 
tion to the ones collected in the real vector x of the primal pro- 
gram. Hence, the effect of the additional constraint 1^2} in 
the dual program on the primal program is, that the set of sep- 
arable states is effectively enlarged. This was to be expected, 
since the set of witnesses is further restricted, and hence less 
entangled states can be detected with further constraints. 

If the dual program finds a restricted witness Z^ such that 
Tr[Zr7] < 1 for the covariance 7 in question, then the ex- 
perimental proof of entanglement can be simplified, otherwise 
some or all of the additional constraints have to be dropped. In 
the program, for each constraint a matrix A has to be defined 
of the dimensions of the covariance matrix. As already men- 
tioned, if no constraint is desired, then constraints = 0. 
Otherwise, 

constraints= [Al A2 ... Ak], (104) 
when k extra constraints of the form of Eq. ( I102> are included. 



IX. SUMMARY 

In this paper, we have provided a complete picture of all 
tests for continuous-variable entanglement which are linear 
in variances of canonical coordinates as solutions of certain 
semi-definite problems. This framework has turned out to be 
applicable both in the bi-partite case, as well as for the various 
separability classes in the multi-partite setting. We moreover 
classified all product criteria, leading to curved witnesses, so 
non-linear witnesses, curved towards the set of separable co- 
variance matrices. 

Finally, we presented the functioning of the two routines 
FullyWit and MultiWit, which deliver just such optimal 
tests, given an assumption on how the state should roughly be 
like. We discussed several examples in detail. It is the hope 
that this picture, and also the freely available routines, provide 
useful practical tools in assessing entanglement in continuous- 
variable systems, also in the experimental context. 
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Appendix: Classification of all entanglement witnesses 

Unlike the case of entanglement witnesses based on states, 
all separating hyperplanes of the set of fully separable covari- 
ance matrices can be very clearly characterized for a general 
n-mode system. This is presented in Ref. 1 1 1], and we briefly 
state this result here without proof for completeness only: 

Theorem 1(1^1) Z is an entanglement witness based on 
second moments in the sense of \2(M and \21\ if and only if 



(i) Z > 0, 

n ^ 

(m) ^str[Zfc]>-, 



k=l 



{Hi) str[Z] < - 



(105) 
(106) 

(107) 



holds, where Zk is the block on the diagonal of Z acting on 
system labeled k. 



This characterizes simply all linear entanglement witnesses 
based on second moments. Further, the set of bi-separable co- 
variance matrices is also convex and closed, since if a state 
contains an (arbitrary small) multi-partite entangled part in 
every decomposition, then all the states in its neighborhood 
will also contain such a part. It follows that entanglement 
witnesses can be constructed that detect only genuine multi- 
petite entangled states. The conditions for such witnesses are 



{i) z > 0, 



(108) 



A/ 



(z*) f]str[z(f,)] > i foralU, (109) 

(110) 



{Hi) sii[Z] < -, 



where 7r(fc) is a partition of the n modes into M < n parties 
as above, and Z^J^^^s^ is the block on the diagonal of Z of the 
j-th party of the partition 7r(/c). For excluding bi-separable 
states it is sufficient to consider partitions into just M ^ 2 
parties. 
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